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FOREWORD 


This work was sponsored under Cross Service Order No. H76789 by the 
Aerospace Environment Division, Aero-As trodynamics Laboratory, Marshall 
Space Flight Center, because the National Aeronautics and Space Administra- 
tion maintains installations and conducts activities along the Atlantic 
and the Gulf of Mexico coastal regions --regions affected by tropical 
cyclones . 

The size and complexity of many space vehicles make rapid movement 
impossible and demand lengthy on-pad checkout procedures. Thus, the 
vehicle and much ground support equipment must be maintained in a 
storm- vulnerable configuration for perhaps 30 days before launch. 

Since this study should also find wide application in a number of 
meteorologicaL^organizations , it is being distributed to several offices 
in the National Weather Service, the Air Weather Service of USAF, and the 
Navy Weather Service. 
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Introduction 


The ‘bivariate normal distribution has been used previously in the study 
of tropical storms and/or hurricanes by Haggard and others ( 1965 ) , 

Haggard and others (196T)» and Hope and Neumann ( 1968 , 19^9 5 1970). 

This distribution is discussed in numerous texts and technical papers 
relating to statistics. Appendix I provides references to some of the 
developmental work and reviews the theoretical basis for the distribu- 
tion. Results of tests described in Appendix II indicate that distri- 
butions of tropical cyclone movement vectors when selectively stratified 
can be described by the bivariate normal model. 

The purpose of this report is to summarize some of the information 
contained in observations of tropical cyclones to provide guides for 
forecasters and the many private and government organizations which are 
affected by these storms. The results presented in Appendix III are 
statistical climatologies of tropical cyclone movements stratified accord- 
ing to season (June-July; August; September; October; November-May) , 
geographical location (five-degree latitude by five-degree longitude 
"squares"), and selected time intervals (12-, 24-, 36 -, U8-, 72 - and 
96 -hours ) . 

Copies of bivariate probability tables and applications by Groenewoud 
and others ( 1967 ) are being distributed with this report. These, along 
with the statistical climatologies indicated above, allow the user to 
make probability statements concerning future storm movements for planning 



or decision making purposes. Appendix IV provides the Fortran IV - 
IBM 360/65 program which easily was adapted for use on the CDC 66 OO 
and RCA Spectra 70/^5 computers to provide tropical cyclone strike 
probabilities which will appear in a second paper. The procedures 
and material presented here should not replace present forecasting 
techniques but should be used as a source of additional information. 

Data Source 

The statistics presented here are based on data taken from the HOAA, 

EDS, National Climatic Center's Card Deck 993 (Tropical Cyclone Deck). 

The preparation of this deck was funded by the Commander, Naval Weather 
Service Command, Washington, D. C. The data are, for the most part, 
taken from the charts of North Atlantic Tropical Cyclones presented by 
Cry and others (1959) and Cry ( 1965 ). A complete description of this 
deck is available in a reference manual available at the National 
Climatic Center. The period of record used here is 1899-1969. This 
deck contains the latitude and longitude positions (in degrees to tenths) 
of storm centers at OOZ and 12Z. All movement vectors were calculated 
using the positions at these times. Only storms classified as a "tropi- 
cal storm" or "hurricane" (winds ^ 3^ knots) and originating in the 
North Atlantic Ocean were treated. These will be referred to as "tropical 
cyclones." Movements for the periods when these storms were classified 
as "tropical depressions" or were extratropical are not included. 

Computation of Statistics 
(a) Stratifications 

The data were stratified according to time and location of occurrence. 
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The year vas divided into five seasons: Jime-July; August; September; 

October; and November-May . This classification separates periods which 
tend to exhibit different characteristics in storm movement or in the 
geographical location of storm development. 'Geographical stratification 
was achieved by dividing the North Atlantic and adjacent areas into 
separate five-degree latitude by five-degree longitude area& or "squares." 

Figure 1 shows these squares and illustrates the scheme used to identify 
them. The three or four digit number plotted in each square gives the 
coordinates of the southwest corner of the square. The last two digits, 
when multiplied by five , give the longitude in degrees . The preceding 
digits give the latitude in degrees. For example, the four digit nimiber 
1010 indicates the area between 10 and 15 degrees north latitude and i+5 
and 50 degrees west longitude. That is, 10°N and 50°W locates the south- 
west corner of the square. 

(b) Coordinate Transformations 

The latitude-longitude positions of the storm centers were transformed 
into positions in the orthogonal T,J grid system currently used at the 
NOAA, NWS, National Meteorological Center. This grid consists of a 
square grid superimposed on a polar stereographic projection of the 
Northern Hemisphere. The transformation equations are: 

I = B[sin (A*)] + 2h 
J = Bicos (A' )J + 26 

where 

A' = (A + 100)(tt/180) 

<!)' = ((J)) (ir/lSO) 
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X = longitude (degrees) 

(() = latitude (degrees) 

B = 31•20i^3 [cos (>p')/(l + sin (((>'))] 

This grid eliminates the curvature effects present in a latitude- 
longitude system. ■ 

Figure 2 shows the I,J grid overlaid on a polar stereographic map of 
the North Atlantic Ocean and surrounding area. The following table 
gives the approximate distance equivalent to one grid length for 
various latitudes. 

Latitude One Grid Length (Kilometers) 


10°N 

2kl 

20°N, 

2T3 

30°N 

308 

i|-0°N 

33I+ 


Conversion from grid intervals to kilometers . 

( c ) Computations 

Movement vectors in terms of (l,J) coordinates were compiled for 
elapsed times of 12, 2k, 36, k8, J2 and 96 hours. All the movement 
vectors originating in a given s'quare were translated such that the 
initial positions were moved to the center of the square. The bivariate 
statistics were computed for the stratifications indicated previously 
by utilizing the machine program. Winds Aloft Summary (1963). 
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The pertinent statistics, including both polar and component forms of 
the means, are listed in Appendix III. These are: 

(1) Resultant direction of storm movement (degrees) ~ (6) 

(2) Magnitude of the resultant storm movement - (D ) 

r 

(3) and (U) Means of the components of storm movement - (AI and A J ) 

(5) and (6) Standard deviations along the major and minor axes of 

the distribution - (s and s. ). These are called SIGX and SIGY 

a b 

in the tables and example applications by Groenewoud, Hoaglin, 
Vitalis and Crutcher (op. cit.). 

(7) The angle of rotation measured counterclockwise from the 
I axis ~ (ijj) 

(8) The number of observations - (n) 


These parameters were computed from the following expressions; 

n n 

0 = Arctan E AI./ E AJ. 

i=l ^ i=l ^ 


= /[( E AI. ) + ( E 


1=1 


i=l 


2 2 
AJ^) ]/n 


n 

AI = ( E AI.)/n 
i=l ^ 


n 

A J = ( E A J . ) /n 
1 

1=1 
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b 2 


^ = 1/2 Arctan ] 


vhere 


(a) AI = I “ and AJ = - J where the subscripts o and f 

or r o 

indicate the initial and final positions, respectively. 

(Note the reversal of I and in the formulation of AI. 

o f 

This modification makes the signs of the components agree 
with the standard meteorological coordinate system.) 


(b) 


^AI 


(( E AI^)/n-l) - (( E Al)Vn(n-l)) 


i=l 


i=l 


is the standard deviation along the I axis) 


(c) 


AJ 


(( E AJ^)/n-l) - (( E Aj) /n(n-l)) 
i=l i=l 


(s^j is the standard deviation along the J axis) 


(d) 


'AIAJ 


n n n 

(n( E AI.AJ. )) - ( E AI. )( E AJ. ) 
i=l ^ ^ i=l ^ i=l ^ 


/n(n 


^^AIAJ correlation coefficient of the I and J components) 


(e) and , the eigenvalues, are the roots of the determinant 


®AI ^ 


^AI^AJ^AIAJ 


with being the larger. 


^AI^Aj’^AIAJ 


s.T - K 
AJ 
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The determinant expanded is 


(s,^^ - K)(s,^^ - K) - = 0 

AI AJ AI AJ AIAJ 


or 


K = 


( s ^ + s ^ ) + y[s ^ + S - 4s ^ ( 1-r I /2 

^ AI AJ ^ ^ AI AJ ^ AI AJ AIAJ 


Examples 

Appendix III contains a listing of the bivariate statistics needed to 
define the distribution of storm movements. Here^ each page contains 
the movement statistics for tvo squares. All seasons and time intervals 
are included except for cases -with less than five observations. 

Figure 3 illustrates how the statistics are used to construct proba- 
bility ellipses. Here, Square 2518 (the north central Gulf area) is 
considered. The data show the end point of the 2J4-hour movements when 
all originate at the center of the square. The season is September. 

The statistics computed from these data are: - (see page III-27) 

n - T3 

H = - . 20 

AJ = -.93 

s =1.08 
a 

% = .65 

ip = 3.5° 
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Longi tude ( ° W ) 



10 


Figure 3 'Twenty- four hour movements ( 1899“ 1969) Tor tropical cyclones initially located within 
square 2518. All initial positions were translated to the center of the square. The 
positions after 2^1- hours relative to this "common origin" are indicated. Ellipses 
for the , 25 ^ .50 and .75 probability levels are shown. 




The probability ellipses are constructed through the following steps. 
Since the statistics were computed in the I,J coordinate system, this 
system must be used in the steps indicated. 

(1) Locate the mean of the movements. 

(2) Construct an I,J coordinate system such that the origin 
coincides with the mean. 

( 3 ) Rotate this coordinate system counterclockwise from 
the I~axis through the angle i|;. 

{k) Select the probability value desired. 

( 5 ) Select the appropriate multiplier from Figure h. 

(6) Multiply s and s (the standard deviation along the 

a 0 

major and minor axis) by this multiplier. 

( 7 ) Let these distances define the length of the semi -major 
and semi -minor axis . 

(8) Construct the ellipse described by these distances. 

In Figure 3 the .25s .50 and .75 probability ellipses are drawn. The 
mean movement vector, along with the size, shape and orientation of 
the probability ellipses, gives a clear picture of how the storm move- 
ments are distributed. Here the data indicate a large variation in the 
direction of movement. The probability ellipses relate the same infor- 
mation by the east-west orientation of the major axis. 

Additional Publications and Future Work 

This is the first of a series of publications dealing with tropical 
storm movement statistics and strike probabilities. The proposed 
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Multiplier (M) 

Figure 4 Radii for various probability ellipses. The 
major and minor axes are given by and 

a (M) . For exaii5)le, the multiplier (M) for 
a .50 probability ellipse is l,l 8 . (Adapted 
from National Weather Records Center, Winds 
Aloft Siammary, I963) . 
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titles and sponsoring agencies for these future publications are as 
follows : 

(1) Atlantic Tropical Cyclone Strike Probabilities (For Selected 
Stations and the Month of September) - Aerospace Environment 
Division, Aero-Astrodynamics Laboratory, Marshall Space Flight 
Center, NASA, Huntsville, Alabama. 

(2) Atlantic Tropical Cyclone Strike Probabilities, (Volume I, 

2ij— Hour Movements; Volume II, 48-Hour Movements; Volume III, 
72-Hour Movements) - Commander, Naval Weather Service Command, 
Washington, D. C. 

( 3 ) Atlantic Tropical Cyclone Mean Vector Charts - Commander, 

Naval Weather Service Command, Washington, D. C. 

Future work may be extended to: 

(a) Use of Hotelling's test to delineate areas of similar 

or dissimilar storm movement in time and space. 

(b) Development of a theoretical model to permit use of prior 
conditions . 

(c) Development of classification and die crimination or clustering 
techniques to isolate homogeneous time-space groups. This 
will be an extension to (a) above. 

SiiTmnary 

The bivariate normal distribution is used as a model to describe the 
movements of tropical cyclones for stated periods from specified posi- 
tions (see Appendix II). Due to the small number of cases, the distri- 
butions are described better by the bivariate t-distribution. As the 
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bivariate normal distribution is approximated by the bivariate "t" 
with an increasing niuaber of observations , it is ass-umed that the 
bivariate normal distribution model can be used to compute valid move- 
ment statistics and strike probabilities. 

The bivariate statistics of tropical cyclone movements are computed 
and presented. Sample sizes range from 5 to almost 100. Obviously, 
more confidence should be placed in those statistics which are based 
on the larger sample sizes. Strike probabilities may be computed by 
the user from tables which accompany this paper as a separate publica- 
tion or by means of an electronic computer program included as an 
appendix . 
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APPENDIX I. THE BIVARIATE NORMAL DISTRIBUTION 


Bravais (l8h6) provides the first extension from the univariate to 
bivariate distribution. Maxwell (l859)» Bertrand (l888), Pearson 
(1900), and Strutt (1919) provide further extensions, and Bartlett 
(193H) discusses vector representations in samples. The following 
discussion is taken in part from Crutcher (l959)« 


A vector distribution is said to be normal if the probability density 
has a maximum at some point and falls off in all directions as 

f(x,y) = exp (-%Q) (l) 

where Q is distributed as with v degrees of freedom for the 
v-dimensional distribution. For the 2-dimensional distribution 


Q= [1/(1-P^")] 


[(X-u )Va ^]-[2p (X-y )(Y-y )/a a ]+[Y-y f/a (2) 

X X xy X y x y y y 


and is distributed as x^ with 2 degrees of freedom, x^ “laY "he used to 
replace Q. The probability that a point lies inside the ellipse for a 
specified x^ is then F(x^<X^p) ^ given probability P, x^^ can 

be determined. Then /x^ or x becomes the vector radius to determine 

^ p '^p 

the probability ellipse contour corresponding to probability P. 


Eq. 1 then represents a bivariate normal distribution where v is 2, X 


and Y are orthogonal components , y^ and y^ are the respective means of 
the components, and are the standard deviations of the respective 
components, and p^ is the correlation between the components. 
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If eq_uals 1, the distribution is a degenerate bivariate distribution 
which is not encountered in practice. The opposite extreme occurs when 
the variances and are equal and equals zero. The expression 
(l) reduces to exp where 

= (X-y + (Y-y and 2a ^ = 2a ^ = a ^ 

The distribution is then circular. These two form the limits of the 
distribution 5 that is, the straight line and the circular. Since the 
correlation between components is often zero, the circular form frequently 
will be encountered. 

Now, if w ^ equals (X-y )^/o ^ and w ^ equals (Y-y )^/a expression (l) 

X! y y ^ 

reduces to 

f(x,y) = exp {-[l/2(l-p ^)][w ^-2p w w +w 2]} (3) 

and if p is zero reduces to 

xy 

f(x,y) = exp {-1/2 [w^ 2 + (li) 

Letting Eq.. (li-) ‘becomes 

f(x,y) = exp (-w2 /2) (5) 

which is the familiar central Rayleigh (Strutt) distribution (1919) if 
only distribution of the magnitudes is considered and the vector mean 
is zero. 

If the distribution is elliptical, then p^y may be significantly 
different from zero. In this case the axes may be rotated through the 
angle ij; to a new axes along which the components are not correlated. 
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The values for the components in the new coordinate system may he 
obtained from Equations (6a) and (6b). 


X' = X sin \p + Y cos T6a) 

Y' = Y sin - X cos \p (6b) 

while the means may be expressed as 

X' = X sin + Y cos \jj (6c) 

Y' = Y sin if) - X cos \p (6d) 


Here i|) (measured counterclockwise from the positive X axis) is given 
as 

ij) = ( 1 / 2 ) Arctan [2p a a /(a ^)] 

Standardization of the new variates X' and Y* provides Equation (5) as 
Equation (j) 

f(x,y) = exp (-(w')^/2) (7) 

and is a measure of the standardized resultant of the X* and Y’ 
components. Thus, the mean of a normal vector distribution coincides 
with the point of maximum probability. In standardized form, the 
probability is proportional to exp I-(w^/2)]. 

Expression (l) is completely defined by five parameters: two means 

(y and y ), the two variances (a ^ and a ^), and the correlation 
coefficient (p^). Moreover, these parameters define the probability 
density as a function only of the vector variable. 
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APPENDIX II. DETERMINATION OF MODEL FIT 


A. Determination of Fit to the Bivariate Normal Distribution . 

This section describes the testing of the validity of the assumption 
that tropical cyclone movement distributions are bivariate normal. 
Crutcher (195T» 1958) made this assumption in work on extra-tropical 
cyclones. Here the assumption was supported by the demonstration that 
the component distributions in themselves were distributed normally. 
Though this is a necessary condition, i.e., that the marginal distri- 
butions be distributed normally, it is not a sufficient condition. It 
may be inferred from Hald's (1952) suggestion (page 602) that a two- 
dimensional test may be made. This was the ba,sis for the assumption 
of bivariate normality for wind distributions as used by Crutcher (1959). 
The reasonableness of this assumption is evident when the expected 
frequencies are compared with observed frequencies. 


Though it may be advisable at times to go to the uncorrelated forms for 
puipjoses of this test, the general case in which the correlation is not 
zero may be used. This is Q or obtained from expression (l) in 

Appendix I. It is repeated here. 


.2 _ 


= tl/(l-r 2)3 


xy 


[(X-X)2/s^2]_[2r^(X-X)(Y-Y)/s^s^] + [(Y-Y)2/Sy2] 


( 1 ) 


where the sample estimates of the parameters replace the population 
parameters . 


Now, the use of the normal distribution implies that a relatively large 
number of observations was available. This is not always the case in 



tropical cyclone data stratified ty season and by five-degree latitude 
by five-degree longitude sq^uares. Therefore, the bivariate t-distribu- 
tion model was investigated also. 


B. Determination of Fit to the Bivariate Student t-Distribution . 

The rationale here is that if the tropical cyclone movements are 
bivariate t and as the bivariate t asjnnpt otic ally approaches the bi- 
variate normal, the non-rejection of the t-distribution would permit 
the assumption of bivariate normality in the computation of storm strike 
or target strike probabilities. The multivariate t-distribution also 
approaches the multivariate normal distribution asymptotically just as 
in the univariate and the bivariate cases. The multivariate form is 
indicated for the ^^-distribution by Krishnaiah and others (1969), 

Steffens (1968), John (1961), and others. Let , x^ , . . . ,x^ be 

distributed jointly as a v-variate normal with zero means, common unknown 
variance , and known correlation matrix ^ = (p. .). Let vs^/o^ be a 
chi-square variate with v degrees of freedom distributed independently 
of x^ , x^, x^,...,x^. Then the joint distribution of t^, t^, t^,...,t^ 

where t. = x./s is known to be a central v-variate t-distribution, Dunn 
1 1 

and Massey (1965)- 


Let random variables x,y have a bivariate normal distribution with means 
and variances respectively, then vs^^/a^^ and vs^^/o^^ 

both are independent of x,y and have a distribution with v degrees of 
freedom where s ^ and s ^ are estimates of and respectively. It 

follows that t^ = (x^-y^)/s^, where x^ replaces y^ and i = 1,2 and each 
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has a Student t-distrihution. The joint density function following 
Steffens (1968) is 


f(ti,t2) = (1/2 ti) 


-(v+ 2)/2 


1 + (t^2/^) + (t^^/v) 


( 2 ) 


Probabilities associated with this function may be evaluated for v 
degrees of freedom and various values of t using the tables developed 
by Steffens (op. cit.). Critical values of t also have been tabulated 
by degrees of freedom. Values of t for a given probability level are 
determined by interpolation using the values of Steffens' Integral 
and his tabular data. The expression 

= (1 - P)A (3) 

where P = probability level, gives the proper value to use in determining 
t when values of have been plotted against t. For example, using a 
probability of .ko and 75 degrees of freedom 

= (1 - .4o)A = .15 (4) 

and interpolation in Steffens' tables yields a value for “ ^2 ~ *904. 

C. Testing of Models for Tropical Cyclone Movement 

Figure II-l shows ten geographic five-degree latitude by five-degree 

longitude squares in the southern North Atlantic and Gulf of Mexico areas. 

These areas were selected to test the bivariate normal and ^-distribution 

function models for the 12-hour tropical cyclone movements during September 

for the period 1899-1969 • The selected geographic areas are shown in 

black. 
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Figure II — 1 Tests of null hypotheses for the bivariate nonnal and bivariate t- 
distribution for 12-hour movements of tropical cyclones during the 
month of Septeiriber. The geographic areas are shown in black. The 
number of tropical cyclone movements is shovn in the upper left 
corner of each square. A single asterisk or double asterisk indi- 
cates rejection of the null hypothesis for the bivariate normal and 
the bivariate t- distribution respectively. The rejection level of 
a = 0.05 involved k degrees of freedom as ten equiprobability inter- 
vals -were selected as class intervals. Period of Record 1899*1969* 
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A standard model was used for testing all squares, i.e., random variables 
AIjAJ were standardized resiilting in means 0 and variances 1. A rotation 
of axis was perfonned to remove correlation. By definition of the t 
variate, these standardized random variables have a joint distribution 
which is the bivariate t-distribution. 

A test of goodness of fit was made for each of the bivariate normal 
and bivariate ^-distributions . For more details the reader is referred 
to Cmitcher and Falls (l9Tl)' The general procedure is the following. 

The distribution is set up with ten shells, each shell holding, theoreti- 
cally, ten percent of the volume. The shells may be rectangular, square, 
elliptical or circular. Availability of polar tables for the normal 
distribution and the availability of rectangular tables for the t-distri- 
bution permits the use of elliptical cylindrical shells for the first 
and square cylindrical shells for the second. The expected frequencies 
for each shell then are n/10 and may be expressed as E^, . An actual count 
of the end points of the observed vectors falling inside each shell then 
is made. This may be expressed as 0^. The difference, (O^ - E^), is 
squared and the square is divided by E^. This is done for each shell and 
the ten quotients are added. This is expressed as 

10 

X2 = S 

i=l 

The quantity is distributed as x^» Pearson (l900). The bivariate 
frequency surface is fitted with two means, two variances, one correlation, 
and a fixed volume, causing a loss of six degrees of freedom. As there 


(0. - E. )2/e, 

X 11 


(5) 
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are ten shells and six degrees of freedom are lost, is distributed 
as with four degrees of freedom. 

Figure II-2 shows eguiprobability ellipses and rectangles of 0.i^0 and 
0.50 for the normal and t-distribution , respectively. The September 
12-hour cyclone movements are indicated by the dots from the intersection 
of the I,J coordinates at the center of the Square 2512. With 30 tropical 
cyclone movements and ten shells, three dots are expected in each shell. 
There are two in the elliptical shell and two in the rectangular shell 
where boundaries are O.ilO and 0.50 probability rectangles. The contribu- 
tion of each shell to for each distribution is (3-2)^/3 or 0.333. 

This is done for all ten elliptical shells or rectangular shells, then 
the total X^ is found for each case. This then is compared against the 
appropriate decision criteria for a = 0.05 with four degrees of freedom. 


The null hypothesis 5. tested against the alternate 

hypothesis H : x^ > \ where a = O. 05 . Here x^/ s is 9-^88. When 

a (a, 4 ) (a, 4 ) 

the x^ statistic obtained is less than the null hypothesis that 

the bivariate normal distribution shows a reasonable fit to the actual 
data distribution is not rejected. 


Table II-l gives the results of the testing of the null hypothesis for 
the two distributions. An asterisk denotes rejection of the null hypo- 
thesis for the bivariate normal distribution, while a double asterisk 
indicates rejection of the null hypothesis for the bivariate t-distribu- 
tion. The bivariate normal distribution model is rejected five times, 
while the bivariate t-distribution is rejected twice out of ten. The 
asterisks are shown also on Figure II-l. 
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Figlire II— 2 Distribution of tropicsal cyclone 12-hr movements for September, 
years l899"“1969 in the I, J grid system. The and .50 proba- 
bility ellipses and rectangles for the bivariate normal and the 
bivariate Student t^distribution respectively are shown. The 
probability of a tropical cyclone occurring vithin the bands de- 
fined by the ellipses or the rectangles is .10. The number of 
movements is 30 . 


II-T 




Both models are rejected in the Square 2018 Just north of Yucatan and 
between Yucatan and Cuba. Examination of the data indicates some bi- 
modality which is evident when the five-degree square is broken down 
into two and one-half degree squares. Tropical cyclones, if in the 
northern part, tend to move north, while those in the southern part 
tend to move west. 

In Square 251^ the slow moving storm of September 18-21, 196k, con- 
tributed greatly to due to several movements in the 0.30 to O.UO 
probability bsind. Though the following is conjecture, this could be 
due to estimates of movement being equalized by the analyst over 
several periods. 


Table II-l. 

Chi-square 

(x^) test for fit 

of tropical cyclone 

12-hour 

movements during September. Period 1899-1969, a = 0.05 

degrees 

of 

freedom = 

critical value of 





y} 




No. of 

Bivariate normal 

Bivariate "t" 

Square 


Obs. 

(rejected*) 

(rejected**) 

1512 


69 

lh .91 * 

5.93 

2018 


5k 

15.60 * 

Ik. 52 ** 

2512 


30 

8.00 

9.33 

251k 


T6 

6.37 

16.89 ** 

2515 


75 

11.00 * 

7.27 

2518 


76 

3.1+7 

8.21 

3012 


31 

7.38 

8.68 

30lit 


Ik 

11.68 * 

1+.92 

3512 


37 

13.1+3 * 

8 . 1 I+ 

351k 


36 

5.67 

6.22 
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Table II-2 gives the results of testing the same null hypothesis for 
selected samples for time intervals greater than twelve hours. Here 
data from all seasons except November-May are used. The bivariate 
normal model is rejected in three of the four cases, while the bivariate 
"t" model is rejected only once. In general, these results agree with 
those indicated by Table II-l. 

Table -II-2. Chi-square (x^) test for fit of selected samples 
of tropical cyclone movements for time intervals greater than 
12 hours. Period 1899-1969, a = 0.05, degrees of freedom = U, 
critical value of ^ = 9.^88. 

(a, 4) 




Square 

Time 
(Hrs ) 

Season 

No. of 
Obs . 

Bivariate Normal 
(rejected*) 

Bivariate "t" 
(rejected**) 

2518 

2 h 

June -July 

37 

5.43 

4.35 

3015 

36 

September 

54 

15.99 * 

8.96 

1516 

48 

October 

44 

15.09 * 

21.91 ** 

2015 

72 

August 

39 

11.51 * 

5.36 

Table 

II-3 provides the approximate ; 

probabilities that the 

computed x^ 

values 

given in 

Tables II-l 

and II-2 

would be exceeded by 

chance. Here , 


for example, the probability level <.05 indicates less than .05 but 
greater than or equal to .02, while <.50 indicates less than .50 but 
greater than or equal to .kO. 


Inspection of Table II-3 shows that if the null hypothesis is tested at 
the .10 level of significance, the number of squares accepted by the 
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Table II-3. Probability levels indicating the likelihood 
that the computed values provided in Tables II-l and 
II-2 would be exceeded by chance. 




Square 

Time 

(Hrs) 

Season 

No . of 
Obs . 

Bivariate 

Normal 

Prob. 

Level 

Bivariate 

"t" 

Prob. 

Level 

1512 

12 

September 

69 

14.91 

1 — 1 

0 

V 

5.93 

0 

00 

V 

2018 

12 

September 

54 

15.60 

<.01 

14.52 

<.01 

2512 

12 

September 

30 

8.00 

<.10 

9.33 

0 
1 — 1 

V 

25lh 

12 

September 

76 

6.37 

<.20 

16.89 

<.01 

2515 

12 

September 

75 

11.00 

LT\ 

0 

V 

7.27 

<.20 

2518 

12 

September 

76 

3.47 

0 

\ T \ 

V 

8.21 

A 

H 

0 

3012 

12 

September 

31 

7.38 

A 

ro 

0 

8.68 

A 

H 

0 

3014 

12 

September 

74 

11.68 

A 

0 

4.92 

<.30 

3512 

12 

September 

37 

13.43 

A 

0 

H 

8.14 

0 

1 — 1 
V 

3514 

12 

September 

36 

5.67 

<.30 

6.22 

A 

ro 

0 

2518 

24 

June- July 

37 

5.43 

A 

U) 

0 

4.35 

0 

V 

3015 

36 

September 

54. 

15.99 

A 

0 

H 

8.96 

A 

H 

0 

1516 

48 

October 

44 

15.09 

<.01 

21.91 

A 

0 

H 

2015 

T 2 

August 

39 

11.51 

OJ 

0 

V 

5.36 

0 

on 

V 
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bivariate normal model would be comparable to the number accepted by 
the bivariate "t" model. In general, however, the bivariate t-distribu- 
tion provides a better fit to the, data. As the bivariate t-distribution 
is asymptotic to the bivariate normal distribution, the difference be- 
tween the two models can be largely attributed to the limited amoimt of 
data. 

In a few squares, topography and/or certain preferred patterns in the 
general circulation may result in a heterogeneous sample within a five- 
degree square. Further stratification could have eliminated the problem 
of heterogeneity but would reduce the sample size, hence reduce the 
significance of the resulting statistics. 

The information presented here is considered to be substantive that the 
hypothesis is valid though it is admitted that the test is not a 
powerful test. Since the Kolmogorov-Smirnov Test is not applicable to 
the multivariate case, other tests are being devised and will be 
published in the paper being prepared by Crutcher and Falls (op. cit.). 

The assumption then is made that the tropical cyclone movements may be 
described by the bivariate normal distribution. To the extent that 
these assumptions may not be quite valid, the tropical cyclone movement 
statistics and strike probabilities will be in error. However, these 
are expected to be an improved approximation to future storm movements 
over those implied by empirical probabilities. 
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APPENDIX III. BIVARIATE STATISTICS OF NORTH ATLANTIC 


TROPICAL CYCLONE MOVEMENTS (1899-1969), (l,J) COORDINATES 
Explanation 

The "5 DEG ID" indicates the square location (see Figure l). All 
distances are given in units of one grid length in the I,J grid. The 
parameters are identified hy the abbreviations in the left column and 
are define'd in the text (pages T“9)- The signs of the component 
movements are as follows; 


I (+) 

I (-) 

J (-) 

J (-) 

I (+) 

I (-) 

J (+) 

J (+) 




+J 


For stratification with less than 5 observations, the statistics were 
not computed. For these cases, zeroes are listed for all parameters 
except the number of observations. 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I>J) COORDINATES 
LDNa 20- 25W SEASON* SEPTEMBER 5DEC ID* 1006 LAT» 10-15N LON* 25- BOW SEASON* AUGUST 
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bivariate statistics of north Atlantic tropical cyclone movements(Ib 99 -i 969 ) coordinates 

L 0 N» 30 - 35 W season* AUGUST 5 DEG ID* 1000 LAT* 10 - 15 N LON* 35 - 40 W SEASON* AUGUST 


Z> •OlA'iTCNJm 

aO'Ocnor*>or- 

xtn****«GO 
rsi ao GO oj (M t-4 
<o 

O' 


oC'4-r'iri^O'or-'#- 

3 • r-4 >J- ITJ (\i IfV • >r-t 

DO'^O^O'^CM 

CM 'O >0 CM r-H r-t 
CM 

r- 


z> • a* r- o CO CM • —I 

asffP»0'00CMC0O 

x^***«*^ 

CM m r<i .-I o 
CO 
>»• 


a:or>4c<t>-imocM'4- 

3 • o o r- o ^ 

Q.-irciO'tnotom 

CM m CM rH r-« o 
M 3 

m 


OCCMr^O'i-Ct-^’OOOM' 
3 •cotnmooh- • r-i 
QCOr-4COrHM>CMCM 

CM CM »-C r-l o O 

' 4 ' 

CM 


QCr- 0 '>l-.-lCMM 30 M- 
3 •rMCMcneo^ •r-i 
□ m— lOOCHr-IO 
Xfci«****(n 
CM r-4 O O O O 
CM 


1- Q- Q. -5 Z 

»-«*-»aoxxacD 
o Q tj o oc a 
> > 

I— I-i-*-3U|UJlJ_ll- 

^ ^ O Q O O 

3 3 Z Z 

tOI/l<<tDQOZ 

IUU1U1USKI~Z3 

oco:zxcocA<z 


O'fMOiTtioaoaoa) 
• sr-^omo • r-H 
CO GO M> CO m >o 

Q£\Oa*«*«crt 
UJ CM M3 -O O CM rH 
CO 

z 

111 

CL.M3aOcn>-HOOr>CMCO 
LU • ^ ^ ITk M3 C<3 • >H 
l/lCDC<3CMOrH{MO 

m • • • • • lo 

M CM m in .-H CM rH 

z 

a 

Vi 

< 

lUOCJ'CO^-cnsrCMCO 
(/) •OOCMI'-COCO • rH 

c^ M> in o '<r 00 cn 

u> ***** ^ 

CM K3 m rH rH O 


OO O M3 ■4" CD fC> O' CO 
• CD CM .H 00 CM • .H 
Or-M 3 0 'OI^C- 

in • • • • *-a- 

CM CM CM O rH o 


iTMnccir^oo'r^oo 

•O^-OCMrH *rH 

CD o r~ M) r- m O' 

^ • * • ■ • ^ 

CM rH rH O O O 


ococfir^r-r^oop 

• *rH 

M3 O' to CC| f*l CM rH 

-M* • ■ • • • m 

CM O O O O O 


1- Qu O. ->Z 
cx:(nzz<crHi-cn 
•hmOOZZOcD 
O Q O O OC O 

> > 

-IJ o Q a a 

33 Z Z 

U1(0<<C3C3OX 

UlUlUltUI»V~Z3 

o:o:zz<oio<z 


ocincMinM-M-cc^crjo 
3 • rH rH CO CM O • rH 
O CMr- O ^ rHM* M3 

JZ '4‘ • • • • • 

CM O' CD CM rH 
M3 
O' 


CC:-4-M>r^M3rHO'r-iO 
3 •CCNOO'lTkrH • rH 

aix 3 M>mo'oo^- 0 ' 

X(<3*****'H 
(Mr- M3 rC3 rH o 
CM 


ac 4 -cMooOMhoocMO 

3 • m O M> m rH • rH 

a-J'CMCMOCMM’Ca 

XC*3*****CM 

(Min M- cn rH o 
to 
M- 


a:mM>rHO'M3rHtno 

3 • CM O' m O' « rH 

OCMO'rHCnOOCMrH 

X(^*****f<) 

CM m m CM o o 

M3 

CO 


i'-or-ioM)^(o*o 

• O rH M 3 O' •rH 
O CO >0 *3 O' O to 

□£in*****o 

UJ(MM>M>rHCM(MrH 

03 

Z 

uu 

h- 

Q.O'i^eocMtnrHcoco 

ui • to to o to • rH 

tococMO'O'rHmo 

St* • • • • • o 

• CMin'S-rHfMrHrH 

z 

a 

CO 

< 

uiO'OinrMcoiftoco 

{/I • rH I— I >!■ CO l'~ *rH 

CM M> CM M 3 M* CD O 

■^ * • • • * CM 

CM CO (O rH rH o rH 


incOrHr^m'O- 4 -co 

■ rH M) fO Sj- O' • rH 

o r~ CO (O o ^ 

^ • • • • • CM 

CM CM CM rH rH O rH 


0 :cD(Ml( 3 ^-rHCOmO CM(^•HU^M 3 t 0 ^»(O 

3 •OO'OtOO *rH • <4- r- ^ CM *rH 

Dor-or-incM<t- to to m O' m 3 co cm 

Z X CO • ■ • • • (O CO * • • * •s^ 

m CM CM CM rH o O CM rH rH O O O rH 

rH ,«• 

I CM 

O 


oe:m( 0 (M(MM> 0 'M 30 o'rHinsrcMor^co 

M3 *'4‘(Mh“l'-0 • rH • rH >S- CM rH M> •rH 

I- Q O' CO o to CM rH *4- O' m CO rH m 

^XcM ****«CO CO*****M3 

3 CMrHrHOOO (MOOOOOrH 

CM 


Q 

O I- Q. CV. -> Z 

rH CO Z Z < »H H- ID 

irH *H a a z z a m 

H o a o o Dc: a 

Q > > 

►H l_l_r-.-3UiUJll.IJ. 

3 -J O Q o a 

O 33 Z Z 

lU (0i0<<CC3OOZ 
Q UiUlUimi— h-Z3 

in o:a:zzcoto<tz 


1- CL Cl. Z 
QCCOZZ<t»H|_cO 
•-■•HQoxzam 
o o o o oca 
> > 

l_l_tH-7IUUJtl.ll. 

_l 3 O Q a o 

3 3 Z Z 

(O(n<i<oooz 

ujujuiUii-f-Z3 

CCQCZZlOOXZ 


III -3 



BIVARIATE STATISTICS OF NORTH ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (1,J) COORDINATES 
5DEG ID* 1009 LAT* 10-15N LON* AO- 45W SEASON* AUGUST 5DEG ID* lOlO LAT* 10-15N LDN« 45- SOW SEASON* AUGUST 

IZ HOUR 24 HOUR 36 HOUR 48 HOUR 72 HOUR 96 HOUR 12 HOUR 24 HOUR 36 HOUR 48 HOUR 72 HOUR 96 HOUR 
RESULT DIR 247.3 249,8 252.6 256,0 261,7 266,6 RESULT DiR 246.3 249,2 252.2 255,0 260,9 266.5 
RESULT DIST 1.172 2.359 3.480 4.547 6.789 8.771 RESULT DiST 1.326 2,566 3,744 4,961 7.322 9.443 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) ( I^ J ) COORDINATES 
5DEG IDs 1011 LAT= 10-15N LDN= 50- 55W SEASON* AUGUST SpEG ID* 1012 LAT- iO-lSN LON® 55- 60W SEASONa JUNE-JULY 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899 - 1969 ) (I^J) COORDINATES 
50 EG ID« 1013 LAT* lO-lSN , LON= 60 - 65 W SEASON* JUNE-JULY 5 DEG 1 D= iOlA LAT* 10 - 15 N LON= 65 - TOW SEASON* JUNE-JULY 
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BIVARIATE STATISTICS DF NORTH ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I,J) COORDINATES 
5DEG ID« 1015 LAT« 10-15N LON* 70- 75W SEASON* JUNE-JULY 5DEG ID* 1016 LAT* 10-15N LON* 75- BOW SEASON* SEPTEMBER 
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bivariate statistics of NDRTH ATLANTIC TROPICAL CYCLONE MOVEMENTS (1899-:969 ) (I,J) COORDINATES 
5DEG ID» t017 LAT“ 10-15N LDN» 80- 85W SEASON* JUNE-JULY 5DEG ID* 1506 LAT= 15-20N LON* 25- 30W SEASON* SEPTEMBER 
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BIVARIATE STATISTICS DF NORTH ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I,J) COORDINATES 
5DEG ID» 1507 LAT= 15-20N LONs 30- 35W SEASON* SEPTEMBER 5DEG ID« 1508 LAT= 15-20N LON* 35- AOW SEASON* AUGUST 
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bivariate statistics QF north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899 - 1969 ) (I^J) COORDINATES 
LONa 40 - 45 W SEASON* AUGUST 5 DEG ID» 1510 LAT* 15 - 20 N LDN= 45 - 50 W SEASON* AUGUST 
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bivariate statistics DF north ATLANTIC TROPICAL CYCLONE MOVEMENTS { 1899-1969 ) (I^J) COORDINATES 
5DEG ID* 1519 LAT* 15-20N LON* 90- 95W SEASON* JUNE-JULY 5dEG ID* 2007 LAT= 20-25N LON* 30- 35W SEASON* SEPTEMBER 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969) (I,J) COORDINATES 
5DEG IDs 2009 LAT* 20-25N LONa 40- 45W SEASON* SEPTEMBER 5DEG ID« 2010 LAT» 20-25N LON* 45- 50W SEASON* SEPTEMBER 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I^J) COORDINATES 
5DEG ID= 2011 LAT* 20-25N LDN= 50- 55W SEASON* AUGUST 5DEG ID* 2012 LAT* 20-25N LON* 55- 60W SEASON* AUGUST 
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bivariate statistics of north Atlantic tropical cyclone mdvements(1899-i 969) n^j) coordinates 

5DEG ID» 2013 LAT* 20-25N LQNa 60- 65W SEASON* AUGUST SpEG ID* 2014 LAT* 20’-25N LON® 65- 70W SEASON* JUNE-JULY 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I#J) COORDINATES 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I,J) COORDINATES 
50EG IDs 2017 LAT* 20-25N LONh 80- 85W SEASON* JUNE-JUCY 5DEG ID« 2018 LAT= 20-25N LON* 85- 90W SEASON* JUNE-JULY 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I/J) COORDINATES 
LON* 90- 95W SEASON* JUNE-JULY 5DEG ID= 2020 LAT= 20-25N LDN= 95-iOOW SEASON* JUNE-JULY 
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bivariate statistics DF north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899 - 1909 ) (I,J) COORDINATES 
5 DEG ID= 2509 LAT= 25 - 30 N LPNa 40 - 45 W SEASON* SEPTEMBER 5 DEG ID= 2510 LAT* 25 - 30 N LON* 45 - SOW SEASON* SEPTEMBER 
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bivariate statistics DF north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969) (I,J) COORDINATES 
5DEG ID* 2513 LAT= 25-30N LONs 60- 65W SEASON* AUGUST 5DEG ID= 251A LAT= 25-30N LQN= 65- TOW SEASON* AUGUST 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) {1,4) COORDINATES 
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bivariate statistics qF north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899 - 1969 ) <I^J> COORDINATES 
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bivariate statistics df north Atlantic tropical cyclone mdvement$(1899-i 969) n,j) coordinates 
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BIVARIATE STATISTICS DF NORTH ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899*1969 ) (I,J) COORDINATES 
5DEG ID* 3018 LAT= 30-35N LON« 85- 90W SEASON* JUNE-JULY 5DEG ID* 3019 LAT* 30-35N LON= 90- 95W SEASON* JUNE-JULY 
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BIVARIATE STATISTICS OF NORTH ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I^J) COORDINATES 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) {1,4) COORDINATES 
5DEG ID* 3512 LAT* 35-40N LDN» 55- 60W SEASON* AUGUST 5DEG ID« 3513 LAT= 35-40N LON« 60- 65W SEASON* JUNE-JULY 
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BIV4RIATE STATISTICS OF NORTH ATUANTIC TROPICAL CYCLONE MQVEMENTSa899-l969) (I,J) COORDINATES 
LQN« 65- TOW SEASON* JUNE-JULY 5DEG ID» 3515 LAT= 35-40N LQN« 70- 75W SEASON* JUNE-JULY 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MQVEMENTS(1899-1969) (I^J) COORDINATES 
5DEG IDs 3516 LAT* 35-40N LON* 75- BOW SEASON* AUGUST 5DEG ID® 4010 LAT» 40-45N LON* 45- SOW SEASON* SEPTEMBER 
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bivariate statistics of north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I,J) COORDINATES 
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bivariate statistics DF north ATLANTIC TROPICAL CYCLONE MOVEMENTS ( 1899-1969 ) (I,J) COORDINATES 
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APPENDIX IV 


Program to integrate the bivariate normal distribution over an offset 
circle. 


The following program represents a variable increment 
numerical integration method as applied to the integral 
of an elliptic bivariate normal density over an offset 
circle. It was developed and programmed by Dr. S. Kaufman 
and C. Groenewoud of Cornell Aeronautical Laboratory, Inc. 
It is reproduced here with the permission of the authors. 
Persons wishing to use the program should try to reproduce 
two test cases before applying it to other situations. 
These test cases are; 

Test Case I Test Case II 


SIGX = 

2 

10.0 

SIGY = 

1.6 

2.0 

CH = 

T.O 

5.0 

CK = 

3.2 

7.0 

R = 

5.6)+ 

17.5 

P = 

0.1188)+ 

0.8)^203 
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program bindc 

DIMENSION VQ(l5),G(5/5/l5)^S(4) 

1001 FORMAT (5F10.3) 

1002 format (2X#6H SIGX /2X/6H SIGY ^2X^6H H ^2X/6H K /2X^6H R 
♦ /2X/6H PROB ) 

1004 format (I10,E10,2) 

1005 FORMAT (5F8.3/F8.5) 

1006 FORMAT (3E10.2) 

1007 FORMAT ( 2X^ 5HS ( 1 > =F6 . 4^ 2X^ 5HS ( 2 ) =F6 , 4/ 2X^ 5HS < 3 ) *F6 . 4> 2X^ 5HS ( 4 ) •■F6 , 

«4^2X.»2HPsF6.4 ) 

1008 F0RMAT(2X/5e15.8) 

1009 F0RMAT(2X^ 15HSDMEBDDY GOOFED) 

2001 FORMAT ( 2X/ 2HM = 1 2^ 4X> 2HI » 1 2, 2E 15 . 8 ) 

2002 FORMAT (6E15.5) 

2003 FORMAT (2X/2HW10F6,2) 

2004 FORMAT ( 2X^ 3HI A* 1 2^ 3X^ 3HI B = I 2^ 3X^ 3HF4«E 15 . 8^ 3X/ 5HDE LX*E 15 , 8 ) 

2005 FDRMAT(2X/3HlAsI2^2X^3HIBsI2>2X/5HINDX=I2/2Xi6HINDX2*l2/2X^4HlN0»I 
* 2 / ) 

2006 F0RMAT(2X^5hIPDS»I2^2x>5HSIGX»F4.1>2X/5HSIGY»F4.1/2X^3HCHsF4.1>2X^ 
#3HCK=F4,1//) 

2007 FORMATUX///) 

2008 FORMATIlHl) 

DO 4 KKK=1/11 
DO 3 III=1#3 

3 READ (5/1006) ( G ( II 1/ J J J ^ KKK ) ^ J J J * 1/ 3 ) 

4 continue 
C0NST=10.**8 
RT2=SQRT(2.0) 

CRTPI=0. 3989422804 

READ (5/1004) IPRINT/ERR 

READ (5/1004) NCASE 

DO 700 ICASE =1/NCASE 

READ (5/1001) SIGX/SIGY/CH/CK^R 

WRITE(6/2008) 

DO 10 1=1/4 

S(I)=0 
10 CONTINUE 
R2=R*R 
RR2=R/RT2 
IPDS^l 

14 GO TO ( 15/ 16/ 15/ 16)/ IPOS 

15 bottom = (CK+RR2)/SIGY 
TOP = (CK+R)/SIGY 

GO TO 20 

16 B0TT0M= (CK-R)/SIGY 
TOP = (CK-RR2)/SIGY 

20 SR»SIGX/RR2 
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ELL=SIGX/S1GY 
INO = l 

RADM45=R2-(4.5«SIGY+CK)*«2 

RADM30=R2-{3.0*SIGY+CK)*««*2 

RAD30 =R2-(3.0*SIGY-CK)«*2 

RAD45 =R2-(4.5*SIGY-CK)#*2 

write (6^2006) IPDS/SIGX^SIGY/CH^CK 

VD{l)=(-4.5) 

V0(2)=(-3.0) 

V0(3)= 3.0 
VD(4)= 4.5 
VD(5)= CONST 
M*5 

V5=(CH-RR2)/SIGX 
VA = V5 
GO TO 80 

30 V6=(CH+RR2)/SIGX 
VA=V6 

GO TO 80 

31 IF(RADM45) 32^33^33 

32 M=M+2 

GO TO 40 

33 IF(B0TTDM,GT.'(^4.5) ) GO TO 32 
IFCTDP.LT. (-4.5) ) GO TO 32 
S45 = SQRT(RADM45) 

V7 =(CH+S45)/SIGX 
VA = V7 
GO TO 80 

34 V8=(CH-S455/SIGX 
VA = V8 

GO TO 80 

40 IF (RADM30) 4U42#42 

41 M=M+2 

GO TO 60 

42 IF(BOTTOM.GT. (-3.0) ) GO TO 41 
IF(TOP,LT. (-3,0) ) GO TO 41 
S3C=SQRT(RADM30) 

V9=(CH+S30)/SIGX 

VA = V9 
GO TO 80 

43 V10=(CH-S30)/SIGX 
VA=V10 

GO TO 80 

50 IF(RAD30) 51^52,52 

51 MsM+2 

GO TO 60 

52 IF(B0TTDM.GT.3.0) GO TO 51 
IF(TDP,LT.3.0) GO TO 51 
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T30=SQRT{RAD30) 

V11=(CH+T30)/SIGX 

VA=V11 

GO TO 80 

53 V12=(CH-T30)/SIGX 
VA=V12 
GO TO 80 

60 IF(RADA5> 61^62^62 

61 IF(IPRINT.EQ.I) go to 64 
GO TO 85 

62 IF(BDTTDM.GT.3.0) GO TO 61 
IF(T0P.LT.3,0) GO TO 61 
T45=SQRT{RADA5) 

Vl3s(CH+T45)/SlGX 

VA=V13 
GO TO 80 

63 V14=(CH-T45) /SIGX 
VA=V14 

GO TO 80 

64 WRITE(6/1008) ( VO ( J J ) ^ J J » 1^ 15 ) 

GO TO 85 

80 1 = 1 
MM=M+1 
MMM=M+2 

81 VOV=VO(I) 

IF(VA.LT.VOV) GO TO 82 
IFU.EQ.MM) GO TO 84 

1 = 1 + 1 
GO TO 81 

82 MI=M-I+1 

DO 83 II=1>MI 

VO(MIMM) = VO(MIM) 

83 CONTINUE 
VD( I)=VA 
M = M + 1 

GO TO (84^84^84, 84^84>30>31>34^40/43^50^53/60^63^61)^M 

84 WRITE (6^1009) 

STOP 

85 X=V5 
VB = X 

GO TO 180 

86 CALL ELIPSE ( I POS/ S I GX> S I GYi CH^ CK^ R/ X/ W ) 

87 CALL A8(X/W/IA/IB) 

88 A=G(IAiIB/l) 

B=G(IAjIB^2)*SR +G(IA^IB^3) 

Cl= G( lAi IB/4)*SR +G(IA/IB^5) 
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C =Cl*SR + G(IA^IB/6) 

D1=G(IA^IB>7) *SR + G(IAMB>8) 

D2=d1*SR +G(IA^IB/9) 

D =D2*SR +G(IA^IB/10) 

E= GdA/lB^ll) 

Fl= A«ELL +B 
F2= F1*ELL+C 
F3- F2*ELL+D 
F4= F3#ELL+E 

DELX=( (360. 0*ERR*SR)/F4) *#0,25 
WRITE (6/2004) I A/ I B/ F4/ DELX 

90 X2=X+DELX 
VB=X2 

GO TO 182 

91 IF(INDX.EQ.INDX2) GO TO 92 
INO = 0 

X2 =.VO(INDX) 

92 CALL ELIPSE ( I PQS/ S IGX/ S IGY/ CH/ CK/ R/ X2/ W2 ) 
Xl=(X+X2)/2,0 

CALL ELIPSE ( I POS/ S IGX/ S IGY/ CH/ CK/ R/ XI/ W1 ) 

CALL NOR (W/PHD 

Ys(X*X)/2,0 

QA= PHI *EXP(-Y) 

CALL NOR (Wl/PHI) 

Y=(Xl*Xl)/2,0 

QB=4.0#PHI*EXP(-Y) 

CALL NOR (W2^PHI) 

Y=(X2*X2)/2.0 
QC = PHI #EXP(-Y) 

Q=( (X2wX)/ 6,0)*(QA+QB+QC)#CRTPI 
S(lPDS)s S(IPOS)+Q 
X = X2 
W=W2 

IF(X.GE.V6) GO TO 200 
IF(INO.EQ.O) GO TO 93 
GO TO 90 

93 INDX=INDX+1 
INO = 1 
XDX=X+.001 

CALL ELIPSE ( IPOS/SIGX/SIGY/CH/CK/R/XDX/WDW) 
CALL AB(XDX/WDW/IA/IB) 

GO TO 88 
180 J=1 

1800 VDV=VD(J) 

IF(VB.LT.VOV) 60 TO 181 
IF(J.EQ.15) GO TO 84 
J = J + 1 

GO TO 1800 
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181 INDX^J 
GO TO 86 

182 J=l 

1820 VOV=VO(J) 

IFIVB.lt. VOV) GO TO 183 
IF(J.EQ.15) GO TO 84 
J = J + 1 

GO TO 1820 

183 INDX2=J 
GO TO 91 

200 GO TO (300*400^500^600)^ IPOS 
300 IPDS=IP0S+1 
GO TO 14 
400 IPDS=IPDS+1 
AAA =SIGX 
SIGX=SIGY 
SIGY=AAA 
BBB =CH 
CH = CK 
CK = BBB 
GO TO 14 

500 IPOS = IPOS + 1 
GO TO 14 

600 Z1=(CH+RR2)/SIGX 
Z2=(CH-RR2)/SIGX 
Z3=(CK+RR2)/SIGY 
Z4=(CK-RR2)/SIGY 
CALL NOR (Zl*AAl) 

CALL NOR (Z2*AA2) 

CALL NOR (Z3*AA3) 

CALL NOR (Z4*AA4) 

Pb(AA1-AA2)#(AA3-AA4) 

PRDB= S(l)-S(2)+S(3)-S(4)-P 
WRlTE(6/2007) 

WRITE (6*1007) S(1)*S(2)*S(3)*S(4)*P 
WR1TE(6*2007) 

WR1TE(6*1002) 

WRITE (6*1005) SIGY*SIGX*CK*CH*R*PR0B 

C interchanges in above statement are intentional, see 400-500. 
700 continue 

STOP 

END 


SUBROUTINE ELIPSE (N*SIGX*SIGY*CH*CK*R*X*W) 
RAD2=R*«2- ( X*S I GX-CH ) 

RAD = SQRT(RAd 2) x 

GO TO ( 10*20* 10>20)*N ' 

10 W= (CK+RAD)/SIG'Y; 

GO TO 30 ^ 

20 W= (CK-RAD)/SIGY 

30 return 

END 
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SUBROUTINE AB { Ql^ Q2^ J ) 

Q3«ABSrQU 

Q4=ABS(Q2) 

IF(Q3.GT«4.5) GQ TO U 
IF(Q3,GT«3,0) GO TO 10 
I = l 


GO TO 12 

10 I«2 

GO TO 12 

11 1=3 

12 IF{Q4,GT,4«5) GO TO 14 
IFCQ4eGT.3.0) GO TO 13 
J«1 


RETURN 


13 Js2 
RETURN 

14 J=3 


15 return 
END 


I 




SUBROUTINE NDR(X^PHI) 
IF IX) 20^45,30 
20 Y=ABS(X) 

1=0 

GQ TO 40 
30 1 = 1 
Y = X 

40 IF(Y.GT.lO.O) GO TO 42 
Bl= 0,319381530 
B2=(-0.356563782) 

B3= 1.781477937 
B4=(»1, 821255978) 

85= 1,330274429 
P = .2316419 
T=l,0/(1,0+P*Y) 

Al= B5«T b 4 
A2= A1«T + B3 
A3= A2*T + B2 
A4= A3#T + Bl 
A5= A4*T 
X2=(Y*Y)/2.0 
C = ,3989422804 
Z = C # EXP(”X2) 

GO TO 44 
42 Z = 0, 

44 IF (I.EQ.l) GO TO 50 
PHI=Z#A5 

return 

45 PHI=0,5 
RETURN 

50 PH1=1,0»Z*A5 

60 return 

END 
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c 


c 


PERMANENT DATA 

2.2E-01 

3.3E-02 

3.0E-03 

4.0E-04 

1.1E~05 

1.7E-06 

ls9E”00 

1.8E-01 

2.4E-02 

2.2E-03 

1,OE-04 

9.0E-06 

4.0E-01 

3.6E-01 

1.5E-02 

1.4E-03 

1.0E”04 

8.5E-06 

3.6E-00 

1.3E-01 

4.5E-02 

1.7E-03 

188E-04 

6.5E-06 

1.6E-00 

5.6E-02 

5.6E-02 

2.1E-03 

3.8E-04 

1.4E-05 

6b0E”01 

2.2E-02 

5.4E-02 

2.0E-03 

5«4E»04 

2.0E-05 

5.8E-00 

7.2E-02 

7.2E-02 

9.0E-04 

2.9E-04 

3 . 6E— 06 

2t4E-00 

3.0E-02 

8.7E-02 

l.lE-03 

l,2E-03 

1.5E-05 

Ib9E-00 

2.4E-02 

II.8E-01 

2.2E-03 

1,8E~03 

2.2E-05 

8.8E-01 

1 « 1 E—02 

1.3E-01 

1.6E-03 

li6E-03 

2.0E-05 

1.2E-00 

1 . 2E—00 

1.4E-01 

1.4E»01 

3. IE-03 

3.1E-03 


INPUT 

AND SPEC 

0 

l.OE-06 

2 


2.000 

1.600 

10.000 

2.000 


CARDS 
4.0E-04 
5,0E-06 
2.0E-08 
1.8E-03 
2.2E-05 
9,0E-08 
3,6E»04 
1 e 4E~05 
9e0E-08 
ls8E-03 
2«3i-05 
9«OE-O0 
3 • 8E-04 
1.4E-05 
9«0E-08 
1,5E~04 
1.3E-05 
1.3E-07 
2.9E-04 
3 • 6E-06 
2.0E-08 
1.2E-04 
4.3E-06 
6®0E-08 
1 • OE-04 
9.0E-06 
9*0E»08 
4.4E-05 
6.5E-06 
8.0E-08 
l,2E-00 
1.4E-01 
3.1E-03 


FICATIDN CARDS 


7.000 3.200 

5.000 7.000 


* U. S, GOVERNMENT PRINTING OFFICE: 1971 —748-000/4011. REGION NO. 4 


5.640 

17.500 
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